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I. INTRODUCTION 


Because of threats from highly Mansuveredls high 
performance aircrafts and the need for increase standotf 
ranges, major improveménes are needed in yuidénce and 
contrel capabilities of missiles. 

The high maneuverability of targats, has led to defense 
missiles capable of devslco higher Lift accelerations and to 
more complex control ilaws, able to improve performance o 
well know laws aS proportional navegaczion. 

In order to accomplish these naw raquirtements with large 


standoff reéenges, BseCDULs On S6SYSo2ns sdS.99 ) air 


eo 
rf 
(Dp 
mo 
ct 
ee 
}2 
J 
ay 


engines has been studied 2nd developad in trecen 


qr 
ney 
iD 
Ny) 
rt 
(n 


The advent of airbreathing engines leads natural +o 2 
Semsoidetation of bank-to-turn missiles in o 
moenangle cf attack of the inlats. 


MerOl Laws, kee 


Q) 
QO 


The necéssity of more complex d 
Prtsale Way, <«O the éepplicaticn cr node=n control and =sti- 
° C 


Meeenon <neOry, since more complzts informat: 


t 


eeaces Of missile and targ2* are nas 


(D 
WV 
W 
Eat 
rj 
kK 


informed by sensors ccmmonly in usé ia 


leads +o the use of a airborne computer. 
meee digital contrcl for application *o terminal guidance 
fea bLank-tc-turn mi 


One  GCN=.1TU0US =WOo dins 


mel lowing form: 


X(t)= A(t) x(4)* Bt) ul(t)+ Eg (1.1) 


11 





where the effect of gravix apeeea=s explicitiy in The chird 
term on the righ hand sid2 of axprassion 1.1. 
After the development of a équivaient dis 


c 
the optimal control problem has been solved,using 


om 


a 
Peat ts, equation due +0 the existence of the thir <3rm 
representing to the gravity effect. 

Next, several analysis has bé2n made in order to check 


the egfrect 


whe 


the effects of small and larg¢s roll excursions 


a 


of the sample rate on the systen, and ene effect of the 


[eed pitch angle, ino check the validity ef such 
m 


rc i 
#wo dimernsicnal model, wh aCe Waa sot 


Maes rst. 
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A. INTRODUCTION 


finecehe pessent work the problen sf terminal guidance f 


ty 


° 
long range, bank-to-turn missiles with ramjet engines, usin 


WO 


a digitalized system has b2en investigated. 


The model developed in reference 1 is used as the base 
Bepechn.s work. Maen tis Grdacaliza=.on Of ehac model, an 
Optimal control law was developed. 


Be ASSUAPTIONS 


tt 


Keeping the same asSumptions as in r2i.i,on= has: 
moe Missil is limited to -2g'sS and #15 gts of 
Pemeeeera lOn in the pitch plans,with zero lag.Also its yaw 
meo Dilot has zero tilag, VaNe nej U4 on Slee nta: ns “Zer> 
Sideslip. 

feSssile *hrust exactly cancels iraj. 

Mie engie Gif attacx :S assumed =D be very small,which 
leads to Siem@Gonmteanaaed ASCeLeration  aciing normal =c the 
Maocs=y Vector. 

Meme tasSsile wail not heve £5 roll through a large 
Meee (EF ULtheS ConSiderations wili give cc this at the end 


Smmeene derivation of <=h= control law). 


C. THE CONTINGOUS MODEL 


Using tne same reference frames as in retrt.l, 


ub 


Sbedy =fane With x, axis perallel <c th 
massile axis, Rocha to) weak owl wore emo 1=ht wing ,and 


Pefeeeeev> 2, axis upward. (see fig.2.1) 


ie 





<j 
a 
t—! 
oO 
QO 
4 
Gh 
< 


SeeeuscnreeteakeceWlth XX. axis parailel to the 


WemeOLr,DOSi tive Z, axis polnting Wpwaeds saad) y, axis 
Bemmeing to the left.(see fig. 2.2) 
ioe sg 2.! and 2.2, che angles Dand & are the Eulerian 
roll and pitch angles. 
The state vector is given as 
= 
i ¥, Y, Rey eZ Ze Beez ap | (2501) 


Seeemyes and 2, 2re the componen=s or the Telative target 
Besltion, Y, and Ze are the relative targe* velocity,At*ty and 
Peeeeare the components of target acsceleration,which is #xpo- 
nentially decaying with a time constanz G. 


Ap: 2-4 2.2) 


(Dp 
rw 
5a) 
t 
(t+ 
{j~ 


where D we =] 


a 
Mme COntrol vectcr is given as: 
a = [ac pc | (2.3) 


wheres Ac is the ccmmanded acceleration and Pec is “he 
Sommand2=d roll rat:. 


The nonlinear plant equation is 


Ix 


eeaxX, uw) + G (2.4) 
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or 


YE ¢ 
Aty tAc SIN D O (255) 
-Aty/& O 
xX = Zz. + O 
Atz -Ac cos g EOE E, 
“Atz/G 
Pe 
Meeteg is dgravity's accelarat:on and GO is =he pitch angle 
as seen in f19.2.2 
en=ea=izing and s¢sttiag 
G=2Eq (2.6) 
one has 
O O 
ae sin D O 
_” O O 
e O O (et) 
Atz + Ac (cos $) AP 
oy, -cos p 0 
- Zz ‘oF O O 
| O © 
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7 ° 
B= Jo Uae eo 9- ccs 6 0 0 | (235) 


in eqn. 2.7,w2 have sét 
cos@ cos(B +4P) =cos fh cos4f-sing sin AP 
sin@ = sin(f +4) =sing cosdfptcosg sin ag 


and expanded in 4® which is considered small. 


® , e e A « 
MeowemeassUuM@ing thet Ac,which is actvally the desired 


conzrel Ac, can be expressed in the form of: 


Ac = Aco 1 ff (2.9) 
7 i 
Bee ACO = Ac at t=0, ons has 
mea x + bu+e8 g (2103) 
where 
Co / O O O O O (22 11) 
DO nO OC) OB Aco [7 - tos B 
Coe 7 OO Oo O 
melo O O O ¢f 9 O 
O O O O O if Aco ne sin b 
GenOmtON OC 60 Us O 
OL VOM Se, OCmeO =a O 





ley 


1Q 


eon 
or 
betw 


Case 


ey 


DESCRET E 


MODEL 


fone Oo Al! otc 


cr 


the continuous-time systemn,> 


inte2rrace in 


eon 


3 4 
—_—— 


Sede= =C 


=e 71 scc ret s¢ 


and 


Will b2 considered, 


in 


Sakic os 


Cont aus 
A=eO-)D 
Samplers ana zero-order holders as i 


such case, considering 
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a 


a ales ee 
h2 nas *o 
Sos -ne 


i 


Cenpuesr 


(25.12) 


(22185 


(Qe 


= 


have some kind 


Sommin cat.on 


— ala 


(Ze45) 





Maomacy Wer 6 the Scare of the Syst2m at time £{(k+7} as : 


Cees 
R(t Y= yet) KEV DE BG) am a(t) (2. 16) 
Caco be 
t/ Pee, om) EM) an 9 (t) 
Ce 


where D(t,to) ms che tranSsictlon matrix or he system repre- 
senzed by eqn. 2.13. 
Furthermore, we will considsred that she sampiing 


instants are zqually spaced, or: 


ct 
‘ 


ket ae a Ca) 


a =—ct + T (22-18) 


so one can replace 
ee = kT 


K 


rus , 


K(KT+T)= O(kT+I) x (kT) (ea) 


KT?T 


+ [Beer B(m) dm a(kT) 
kT 
AT#T 
+ |b cerseey ) ECM) dy g( kt) 


KT 


Seeen a Simplified acotation: 
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x(k+1) = Ag (kK) x(k) + By (kK) u (k) 


Seren) og (\) 


Ag (k)= © (kI+T,kI) 


KT+T 
By (k= [ B (xt4T, my) Bln) aq 
kT 
kT +T 
sgn [ Beery ) Bl) am 
KT 


2. Calculat 


-= —_™ 
— —— ow 
= —— =P 


ies Sera =gnerorward te-show, 


Pemee Matrix A(t), that the =ransiction matrix of equation 


2-19 is: 
ee Ad , O 
0) Wie e) 0, 
0 0 e@ O O 

DikT+T) 2 OG 6, aa a 
O CG: 0 a 
OOO. om 0 
Ooo 0 O O 


BeOS. ane eee Suenne): , 


OQ Ad,, 
O Ad,, 
OM 8, 
Ad, Ad,, 
Ad Ad, Y 
eE O 
O { 


3B (kK) 


Iu 
fiz 


jy 


a2) 


G2) 


(23722) 


(22-23) 


(x) 


using the sparseness 


(22-24) 





where: 


= 
oe! = G4 (1-6 @ ) 


Ag . =Ad 


i,3 4,6 


(1-2°% ) 


Ads s aaa yt 


POoterie Calcula=2on, of the others terms one may make 


Meme r cne property of the transictioa matrix that: 
— = A(™,) Aq (t, ec) 


SO, 


Kits 1, 
d Ag (kT, kT) = A(kKT+T) Ag (kT+HT, kT) 


d(kT+T) 


BotcoA (2,7): 


Ady (tet kT) 


1) 


e 


Ay (KT+T) = Aco 1 | Es 


d(kT+ T) tf 
kT + T 
: T: SIS Pp 
d(kT+T) 


KTeT 
i ae KT + T = 
Sd Aco sh, [ | my Oct er) 
KT 
a2K+i 
A = Aco | =e eo 2 
ae Cc cos £ T | 2: T 


nas wate igs) 3 


KT+T, kT 
_dAa,, (kT, k Ol Adz, (kT+T, ey 
d(kt+T) 
«eT +T 
Ad, (KT+T, kT) = [| Ag a (kT+T) 
kT 
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ft Zeal 
Ad, , | | oT | © | Aco cos g 


Dozmnge=n]= Sale process for Ad{5,7) and Aa{(&,7) 


Rom the SSevdie on Ome snow Mack sk oo) (K) 


meeeeaing ©O eqn. 2-22 


OkT+T) B(7)= Ad (KT+T, 4) BCH) 


where, 
/ Am An oF (9, 
Oo f A, 70" eo 
2.3 
O 0 An 0 O 
3,3 
Ag(kt+T,9 ) sn 
Oo MO 8 Oy 
CO Ee (0) 6, 
Oa Omen ae 


where Ay represents Ag (kt+l, ) and 


O 0 

Sin Q, O 

O O 

O 

Bin) : fs g O 
O 0 

O 0 
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o= 


one needs 


O An 
es 
Orace 
aula te 
are Nes 
oe 0 
@ { 





thus 


An Sin A 


Sin Be 


Ag (kT + Tm ) Bem) = “An, <o5 6, 


- CO5 2 


AGO cos B 


a cae " 2 = 
atlay | racren) a (kT+#T) T 
2k-+ i! . 
age [kate] ong 
Gkiee @) 
tes. xT ae. 
al | : ZT i ZTi 


eed «LOT By (k): 


ZZ 




















Ba. Ba (2. 25) 


,2 


5 
i Le Ga,, 
0 0 
By - Bd, Ba, 
- Cos 2 Ba 
O O 
O 1 
where 
kTe T 
tice |S 
Ban ert) dy sin@ = am sin g 
KT 
Té 
a, in cos B 
kT+T 
Bd = A | 
i. | Ver 
KT 
Which after some algebric work has been found 
Port | (ere)? © (eter) at. © Ue 
ee ee aay | Acc cos @ 
Z Soin Za ie 671 2 
and 
kT+T 
an 7 pica 4 
kT 


wmech can be found «o be: 


Z3 





ne i i 


thus 
ees eZ kas | 
Bd, {X) = 3 3 OT: re Aco sing 
2 3 a 3 
$.2 2 3 Ti Zn 6T} = 
In the same way: 
O 
O 
O 
P(kT+T, a ) E(m) = -T cos 8 
-cos § 
O 
O 


meen (K) <s equal to 
Té 
z= | 0 Sono! 0 ry ee -[T 20s 9 0 0 


whezre we have considere? 9 as a constant Se owe (tse ner 


memmen=S Ch this after development of che centrol law). 


Notice that neeugouecmpcioS “Ost ,2h>. econmanded 
acceleration has been considered an unknown and the assump- 


ct 
1° 
O 
iJ 
ao 


G@eepeeci made i<chat it will be in =he forn of eqn.2.9. 
Mies Will nesd in Eurther develspments «9 consider the 
Control Ac as a known Ac(k) ,which will ba a constan= between 
moeand KT+T. With such assumptions, the discrete represen- 


tation of the system is easily round <*> be : 
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x (k+1)= Atk) x(k) + E(k) Ut EG 
where 
Vial Ad , ga0 oO 
O f Ag, Oe Oe 
pu ia OC 
Ag = O O O / if Ad, . 
Canon O. cO Ade, 
-T/, 
Oo @ @ C6 OES 
Guo Oo. Oe O 
with 
Ad, o T- 6° (1- 2%) = Ady. 
ey 
Ad, > 6(1- 2 /e = Ade 
Té 
aie =>- AC cos b 


ae tanec COS 2, 


rT? 
= --- Ac sin 
aes) 2 Pr 


Ad._= 


Bt fen Gc. Sin po. 


ZS 


7 (227) 


Ad 
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Ad 


4.7 


Ad 


x7 





and 


T sing, T he cos 
T sing, = Ac cos B 
O O 
By (k) = - - COS D = Ac sin B, 
aaa cos B = Ac sin B 
O O 
O a 
and 
O 
O 
O 
Eg (k) = -= cos 8 
Simcos: Oo 
O 
O 
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Ac stn D 


~—e «me Re ee —= SER qe qe qe cer aa 


Ac cos D 


Frames. 


Reference 


Figure 2.1 
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me eee a eee eee eee ee ae ee ep ee ee ee SD ee ee ae eer a a Oe ae ee ee 26 Ses See 


ee eS SE es CE eee: ee Se ce ee ee ee ee ee gee SE ee ee ee ee ee Se Ge SS ee Se ee ee 2 Se ee EE SS — 


ne 


Ragune 2.2 


Frames. 


Reference 
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TIL. THE OPTINAL CONTROLLER 


A. DERIVATION OF THE OPTIMAL CONTROLLER 


Moe order to have a suitable guidance law cto inpl 
the control commands,we will Sloe zinGgeeen>s 8 following 
performance index: 


Mat 


J= -5- x"(n)a(n)x(n) + -5- w(k) 2 (Kk) u (k) (3.1) 


kz0 


where x(N)is the final state at t=Ti 
As we want to minimize the final miss distance,the 


weighting matcix W(N) is taken as 


i Ceo mo eo) Oe 
(3.2) 
Geen OC OO} OD O 
Omewoe OFSC OO .0 0 
0 Camo 1 ODO 
O Quam ne Oo O O 
O OO Oo ONO EO ege, 
O Cm iO CO maa 
and Q(x) Foo wOurDy “WO SOc ei Ve “dsftinits synmerric 


Weighting matrix to be chosen. 
meee ne derivation of the soluzion, cxrcefer 
followed keeping in mind that the state squatio 


=orcm: 





xX(k+1)= A(k) x(k) + Bik) u(k) + Bg (3.3) 


x(k+1)= £ (x(k),u(k),g) (3.4) 


Meeee i =he “hird térm,which represents <hs sffect of the 
gravity has b2en considered constant. 


Monsideting that the partctmanc= index is in the forn: 


Nel 
J= D[x(v) ] + | Lx | x(K) ,u (kK) 3 (6) | (3.5) 
k29 


we need to find &@ sequence cf u(k) that minimizes J. 


Micon the system equsticr =o J With a multipiter A (x) 


N-! 


s= fb {xo ) L(x) x(k) ,a(k) 3} + (3.6) 


k:O 


z 
+ Xd (k+1) { 2, [ec 6909 3) : cose | 


and defining a scalar sequ2nce H (k) 


pC kK) = boy | ca 009g | r (3.7) 


r 
+ A (+1) | x (kK) ,u(k) 3 | 


Re Oli ys ose No | 
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one has: 


. 
s= Bf xan) - A(N) xy) + (eo 
Nel 


+ ) [Hod - XK) x08) | + H(0) 


K= 
Mensidezring difilerential changes ino J: 


_| LY | | 
dJ= aa d (N) dx(N) + (3n9) 
Nel 
JH (k) x 
eyo k 
+ ( ] act) + 
k=l 
» VH(k) ducky} + 9H fo) dx (k) + 
du (k) d x(0) 
o u (0) 


Gaeos-ng the nultiplier Atk) SO 1a 


Xk) = o H(k) =0 


Bel 
x(k) 
hus 
dL (k) + en dt Ae ie 
d x(k) 9 x(k) 


sia) 








and 


Bate) _ 


(3. 11) 
od u(k) 


with boundary condition 


A) = 38 (3.12) 
o x(n) 


we obtain the minimization of the per forneance3 


ee Xe, 
In the present case wa have: 
Nol 
] T T pairs 
J= 7 Kernen) < (hn) + US Gk ep Uk) (2-13) 
T KeO 
+ A (k+1) ] a (xy x(k) #3 (KI u({ XK) +E so x(k +1) | 
eeed H (Kk): 
\ 
H(k) = 5- u' (k) 2(k) u(k) + (3. 14) 
he 
+ D (k#1) Ee Cikinet ot eek). eS 2 | 
mimeo an Order to minimize J: 
OH Ge) u? (ky) Q(k) + » (k+ 1) B(k) =9 (3. 15) 
a u(k) 


a2 








= 

or considering that Q(k) =0(k) 
it 

Q(k) u(k) = - B (k) Xd (kt) 


and 


Koy = 2H) 


jag (3.16) 
o x(k) 
T 
\(K) = \ (k +1) A (k) 
and 
T 7 
N(N) = x (NP WN) (S217) 


Mietace that we are not weighting the states in the 
Memeo sNance :21dex, except the last stat3s.A more general tcrn 
ested b= ob*+aine2d, With 200 )Ste82s pS 2aG Ww 


Smenge che ¢qn 3.16 to; 


T T it 
\(kK) = ACK +1) ACK) # x (KY) W1(K) (3.18) 
meets W1(k) is the weightiag matrix of th2 state 


(eee 1S 


Ss 
PeeeeeguaescnS 3.15,3.16 and 3.17 ont 1353 ables t0 find 
meme sequence 5£ u(k) that will give che minima co 

S 


Such set 9f aquatics can be solvad by che sweep nacthod 
Meme) cref.2 . 


Aen OOK £at a Solution of th] Eorn 


ee 


u (kK) =-F (k) x(k) -FG (k)g (kj (3.19) 


aS 





Mea] Means that the commanded acceleration and roil rata, 
will have a correcticn dus to the affect of gravity. 


Placing: 
A(k) = S(k) x(k) + SG(k) g(k) (Be) 


peom eqne 3.15 
r 
QO (k) u (k) =-B (kK) S(Kt1) x(ke 1) #33 (K+ 1) g (X41) | 


aeoml eqn. 3.3 


Q(k) u(k) = -B’ (x) S (k+1) } ack) X(K) r+ (332 1) 


+B (k) u(k) + Z 3 | - Bk) SG(k+1) g(ktl) 
So 
lo cx 03, Ukeyies (K+ 1) 3 | u (K) = 
-B’ (k) S(k+1) Ak) x(k) - B (ky) S(k#1) # g(k)- 


-B’ (k) SG(k+#1) 3(k+1) 


Sensicering that g is a constant 


u (k) =-}o(n + B’ (k) S(k+1) B (%) | (sere) 
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lB (iy S(k+1) A(k) x(K)+ 


+ | Be) S(k+1) E + B’ (k) sa (ket) | 5 | 
SO > 


u (k) =-F (k) x(k) -FG (k) g 


aoe 
[ ocx) + B (k) S(k+1) B(K) | 2 (3.23) 


ry 
it 


- B (k) S (k+1) A(k) 


FG = es t B. (k) S (k+1) slot (3.24) 


{ 3" (k) S(k+1) B+ B'Sg (k+1) | 
mem eqn.3.16 and 3.19 
\(k) =a" (ky) A(K+1) =A (K) [See (ke) #56 (ket) g (ke) 5 


T ky S (e411) ACK) x(k) + a (kK) S(k#1) BK) uk) 


=A 
+a” (k) S(k+1) Bg + A’ (k) SG(K41) g(k+1) 
feom eqn. 3.19 
(kK) = a” (k) S(k+1) A(kK) x(k) + a (Kk) S$ (k+1) Bk). 
» [=F (k) x(k) - FG(k) s | - a (4) Si eee > 
+A SS(k) g (k+1) 
Sen) aS g9 1S a constant: 


S (k) x (k) +SG(k) g= }*" (kK) S(k+1) A (K) -A (kK) S(K41)B(K) (K) | + 
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“\a" Go S(k+#1) E - A (k#1) B(K) FG(k) 


+ A (k) sc (k+1) | g 


+hus 

S(k)=A\ (k)S(k+#1) A(k)-A_ (Kk) S(&#1) B(K) F(K) (3. 25) 
and 

SG (k) =A° (k)S (k#1) B-A’ (k) S (k#1) B(K) FG(k)+* (3. 26) 


+ A (k) SG(k#1) 


Mryese equations, 3.25, 3.26, 3.23 and 3.24 can be solved 
PeeGnwards with the final cond: *ion: 
S (N) =W (N) 


SG (N) =0 


Seq 2foe tos  POUnmy ==VeoNS. DOUG dazy -cond=~ 


(N) =W (N) x (N) 
S (N) x (CN) =W O(N) x (YN) 
SO 


S (N) =W(N) 


Be EXTENTION OF THE MODEL FOR LARGE ROLL ANGLES 


UoecOuecnLS point, See tasecom cake i182o. account *hat 
throughout the development of this work,the angls has been 


Sens. dered small, it is necessary to rslax this restriction. 
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EMmoOnrder tO da thas, the system has been broker in two 
Becks as in figure 3.1. 
fremitirst block 25 2 representation of the elgoritkhn 


which will calculate the optimal connands. The algorithm ha 


™|™ 


4 


contain with itself an exact model of the svstem or nissils. 
imemmodel of the miecsile 1s initialized from the information 
on the initial states, the initial input command Aco and 
[eed ! oll angle (at t=0 ); and computes tha sptimal gains 
pm@eeruccher the optimal commands whica will be fa2ed to the 
missile. 

The method adopted in computing <*h2 cptimal commands is 
more easily understocd if one consiiers figure 3.2. 

In figure 3.2, the lines numbered as 0 in the araph for 
Mememnd £cr Pc are the optimal commands for a given initial 
roll angle(g,). Lines number 1 are th2 commands for a second 
mete cal roli angle(@ ) larger thang, SecewcOmome ois in 
megmce 3.2 one has a family of optimal commands for any 
Merc aal coll angle. 

[fee es =hat =he upper Line of the graph of Ac 
eepomaccelscrations of a missil2 whith had at *=90 ac 
Metal coll angl= in order to hi g 
BemendsS in roli rate. 

Mmieene present method the cemputer pertorns the calcula- 
mom Of “he commands only for tha first st2p of tine and 
then feeds thes2 commands to the nissils. The missil2 is 
then driven to the next state (x(k+1)) and f2ed-ba 
Bemeucer “he information on the roll angle at that step.The 
roll angle fead-backs from the missil2 is considered bv che 
algoritm asthe initial roll angle act t=0 ar h 
Bemmands ars Calculated.Notice that at this sécond step «he 
algorithm will feed to the systam the second command (at 
t=t, ). This process is them repeated until = is equal ‘to 


the intercepted time. 
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Pcp aan ee sOuncal~ze thal with =his methced of 
calculation, Since the algorithm was developed with the 
assumption cf small roll excursions some error is expected 
@ue to the fact that in computing the gains by solving a 
Ricatti equation backwards,as has besn done,it is necessary 
to update the system from t=0 to t=Ti at @ach step, and iin 
this process the roll angle is not small.Notice however that 
w2 are applying the commands only in one s*ep,and if one 
expects that the roll rate@ will decrease to zerc,2zs we are 
[Miezeaqsing in time, the variation of the roll angls will 
tend to decr2ase,sc,we can expectsda that the error will 


decreasé as the time increases. 

Another important point tc be studied is how thes 
mese lt (second biock in fig.3.2 ) has to be implemented in 
Mea== to be valid for small and large roll exc 

It is considered that one has the perfect kn 
the commands, thus the missil2? is a2odéile Sa a Vv 
meee SyStem as in eqn 2.26 with the initiai reil angle being 
update at each step. In this way the system wiil cake <zhe 
Metal roll angle as tha Summation of all previcus initial 
roll angles. (seé fig 3.4) 

BEOMche O=i1dinel variables one has Ecor “the etate 40, 


eee follcwina: 


Ag+) =40(k)+ Bens (gee) 


}- 
t! 
ct 


(D 
"a7 ;) oD 
e : 


iQ 


meen 25 Shcw in rigg. 3.3, where the ial roll angie is 
kepz constant and Ag is update each st 
2 -7oll excursicns 


the 2ng12 @ has 


}-- 
oy 


Notice however that considering 


ct 
a 
jv 
at 


(see fig. 3.4) and keeping in min 


bean defined as 
B= ¢,+49 
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Bie expression 3.27 is not valid, Samee =n  medaling =h- 
systém it was assumed that S® would be samll. 

This leads to a change in the =xpression for the varia- 
ion Of Ag Scere d Gs -654  Etrag. S24, the angle A®is 


update at each step, so one has: 
Ag(k+ 1) = Pc (k) T (3. 28) 


By consideration of the equatlon 2.27 it can b¢ seen 
that by setting the element A(7,7) *t0 zero, one cén obtain 


Peta c.On 3.28. 
Maus the missile modei for large roll excursions will be 


represented by the fcllowing Equation. 


x (k+1)= A (kK) x (kK) + BK) u(k)+ Eg 


where 
/ ne, O O @ A,, 
(32529) 
O ata, OG oO A,, 
-% 

CeO oD) ORO eo O 

A(k) = O O O / a Ais A, ; 
0) OOO ey b. 

% 

1, OO O GQ «°° O 
0 O O O O O Aa, 
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where A(7,7)=0 ,for large roll excursions 


A(7,7)=1 ,for small rell excursions 
ies — Ac cos 
Cee 2c : az 
T?é 
T sin %, aC cos B 
0 O 
(3220) 
3 
i? as d 
T¢ 
-iecos DB — Ac sing 
° Zz 
O 0 
O ae 
O 
rz (33231) 
Ze 
-T cos 9 
0 
O 
Hemera d=erijne the staves *«,, aS the relative posi- 
Peon, X22 aS the reletive velocity in the ¥ direction and x, 
Pen ee aaget aC Cele raclOn in Y ditecciom. The states x, , 
imeand Xx hes the same meaning, ier wecne Z-aiseczeton an 


4 Q 





a 1S ap. 


In «his repres 


equal to the component of 
mepection and: 
pee y= ¥,(ktl)= y. (k)+ 


on 
+ [2 Peeoenti— eG | Aty(k) + 


Z 


3 sing Ac (k) + 


i 


Meer 
conm 
térm represents the effect 


to be small due to the cub2 


Hee the Component of on 
on? has: 
mre iy= Z,(k+ I= z.(k)+ 
+ 


k T-Z2(1- 2 4) 
2 


oe COS @ Ac (kK) + 


entation of the systen, 


X, (kK) is 
the miss distance along “ne Y 
Y, (k) T + (3.32) 
tia : 
eae Ac(k) cos B 4 B(x) 


T3 

as Re CS) cosh Pc(k) 

Geeaew Rho Or SGhw 3.52, ase casiiy 
to the miss distance of respec- 
distancs, relative velocity and 

Ollowing two terms represents the 
and 2deeseselesezaon and the Last 

O-GOUPleodaeec amd PC, “and fends 
Galen Ss aeee “pet io . 

wos ecto e me the 2 07 = Sots On, 

Ze) leas 


2 
Dee Zale -5- Ac (xk) sing A(x) = 


+3 TT? 
aoe entst sing a ic c= Cl S 


whe-¢ its “*érms have tne same physical meaning as in the 
PeaeseSss- On fOr x (k*t1), With the sriect of the gravity added 
to the 2xpression. 

BENGE Oe Can Nnotece that in che representation of the 
Messe dietance in Y direction appeal two cerms as a function 
ef cos @, and in the representaction of the niss distancs 
mies Glrecce=0n appear two terms as function of sin g eee SS 
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interesting to verify that the fourth term in the RHS of 
meme XDLESSioOnsS acts liks a correction for the fLicty tern. 
FeGeting to fig. 3.4b, one can see that at any step of tine, 


the commanded acceleration is actually, 
Ac cos @ -AAc 
and considering small angles: 
A Ac= Ac cos, - Ac cos(g, + Of) = 
Ac cosQ, - Ac | cos sf Ap- sing sin Ap | = 
Ac sing Ag 


The same idea can be applied to the ¢xpression for 
x (k+1). 


emcee ard) X-, | ©Spresens che Telative velocity, 


and are 
@ ade 
macs * 1)= y, (kK) + G (1- 2 Zs Aty(k) + Ac(k) eae T 


aoe 
@eac({k) sin® T+ Ser ciais) cosg Pc(k) 


rT 
% -’%z 


fck+ i= z.(k) + G (1- ¢ ) Atz(k) + Ac(k) sing T - 


2 
- Ac(k) cos B T + 5- Ac (k) sing Percy -er cos 8 G 


Where the two SaSt (SUNS Ln fieekaos.) beprtesents the efrec= 
Bemeene Velocity and acceleration at a previous step, and the 
ether three terme has the samé meaning as pravicsusly stated. 
ea oe cad xX, eateeenewcargse. ‘accelerations, in 
this model being expcnentially decaying. 
1. Ffsczs on the Miss Distance of she Extsnticn oF she 
Model 
In previous subsectior, extention of the dee fos 


large roll ¢xcurssions has be 


een pe m 
are *wo models of the system being used. The first 
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in the algorithm is valid only for samll roll ex 
and a second mcodel, valid for small and largé roli ¢xecure- 
sions used as a representation of the missil:. 

The algorithm with the first model, as 2xplained 


Vole a tizeq 2c “ea 


on 
iD 
th 


4 
id 
- 


€2p with the actual roll 


c S 
enemm 2 SSil, and performs the calculacion of ths 


9 AQ 

» kF O 
ss) OD 

Qu 

in O 

e t h 


Q Pm 
Oo 


Mmorieé= =o check the effect of the extention 
1 


iD 
f 

j y ] 
fu 
(D 
pm 


model on the miss distances, one can iain: 
roll angle( fo 
have the comanded accelaration vactor pointing ¢+ 


O° 
n 
Pees eo ce meet ore ae 2-0 ~ in osdsc to 
° h 
projectedq final target's positicn. PieSeteans =hat the 
Meeete WOUlJ not have to roll tc hit the ctarge= (se 
eases |}, hus the commanded roll rate C E 
fewer ienm Will b2 equal +o zero. This 1 


e 

mpli 

point ahead, the BOLL) “angle gets ‘Constan= and equal 
aces 


fu 


eee. , and that the time history sf the control ac will 


}+ 
® 
6. 


Pema siutaight 1 
mee fact th 
des2rves an inv 


Ra SNOt 8G euthae, pas) Ss Oech Giewes-s.2 5 
m 


W 
Bo Mie ag. ey aS 


ae <~he moment the missi 

Semeerol Ac will be «th Beg “Sed Se@eomess=icn <O hit) the 
Saeoget if “he missile initial roll angle was ee ne ees 
Peremenat ~he contrcl Ac computed would be correct only it 
BpeeMassile would have turned immediatly tec this angle. 

For varv small angles, *he previous comment would be 
acceptable, Diktemeon) a NROXrMaAl situstion, as showed in 
[mee ]ss.5, the missile will only csach the idsal initial 
Boll angl2 after some tims, and jdus 7-90 the vertical ccmpo- 


Meme ot the acceleration, when =his occurs the missile would 
Memes a DOSiticn abcve the ideal trajectory, which m2ans 
Sage i= Would follow a course parallel to the idsaal trajec- 


Beey tO intercep<. 
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Oré can see that such problem will lead +0 a larde 
Miss distance, in the case that th? target accaleration it 
femmon smell. Thus, some corzection is necessary in order <c 
impreve the missile performance. 

Figur2 3.5b shows the missil at some point of its 
trajectory where 1t has reach its maximum roll angle, thus 
pomaes at a parallel course with its ideal trajectory to 
meet cept. hae ies COLA: a stmecuall the states are known, 
it is possibl2 to recompute a new D scat ‘s Sone be ee cols 
Peremre che computsr is fsed with with the states at this 
Domnt 1+ will compute the commands in order tS drive the 
missil to the new Ce PpectecinD tee NtTOducs the desized 
goa rection. Notice that such correction can be made during 


aeeecne flight, 


en 


plished by fesding-back 


all the sta 
computer 
commands at 
being updat 
mie-=rcept, 


intercept, 


where Ti 


*+es Gr 


the present 


zach step 


ed, Ze 0 es 
which has 

fe paged 

=k - 


C. SIMJULATION 


In order 


+ha max~ri3 


Q (weighti 


EL Ole c=O0n: 


dO 


cr) 


+he missile. 


2S able <> perform th-> 


Oe 
also 


been 


1S tne fon nal “2mne 


<0 keep the 


ng 


suggested in such reference 


Oe vata 
eK son as 


the 


nezhod 
compute: 
Wee h 


Caleta On 


> 
SheoSe ne 


time. howeveras ~he states ars 
necess2cy £0 upda=2= the time to 
SOtemUerIg “che <ame co ge Ec 
TOMNac =. GSD cs 
Sam2 assumpzions as reference 1, 
She ween 1) has been put as 


= @ 
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b, 0 (Bes) 


meve different cases were run: 

Cas2 0, *¢sted with one simpie meodei valid only for 
Smelt roll angles begins with missil and target on parallel 
Memm@s]es to the inertial x axis, tha targe= 190 m¢tars above 
the missile and with an e¢vasivS manoever expo 
Beeeeeng with time constant of 23 seconds. The initca 
Meee tat: On Df sargqet was -.5 g's in y direction. (see fig. 
BO ) 

Case 1, with the same scenarios as case 0, but was rua 
Memmog che algorithm fer large rt 

Case 2,is th2= same scanari 


O 
° 

the initial acceleration sf target w2 
= 


a) 
— 
4: 
ey 
4 
ee 


direction. (Same as cas 


ese 3, the Sdm= as caseZ2,with target acceleraticn cf -4 
1 


g's. ( Same as case 2 in tha reference 1) 
Case 4,same as case 3 but with target at initial posi-~ 
tion 600 meters bellow the missile. (Same as cas 3 in refer- 


> 
eénc2 1.) 


D. COMMENTS AND CONCLUSIONS 


1. PBesulzs 
icici se- ths missile bagins its trajectory 
Semmanding 26.5 m/sec@ and the tims history of «he control 


Mm~mpeorelows eXact ly a Straight tine in the form suggested in 
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UJ 
in 


meee 1, a= seen in fig.3.7The control Pc, begining a= . 


rad per second, is decaying and c2achs zero at K=80 (ses 
Paqeeoec). Figures 3.9 and 3.10 shew ths miss distancs, where 
one can see that the missile is crossing the target with a 
CG-to-cG distance of 1.5 meters. 

Figur23. 11, SHowsetne tate history io2 the roll 


angle, which as expect, reachs a constant valus, with the 
missile crossing target at t=T1, with a bank of .44 radians. 
In case 1,which was run with the model for large 


roll excursions, the missile kept “h2 same Aco, but there is 


e 

avery small increase in further commanded accaleration in 
meme -O CcOtrect the effect of the roll angle cn the 
Peencel component of the control Ac, as seen in fig. 3.12, 
and table I. 

Figure 3.13 shows that tas roll rate decreases 
muetest as before, andthe final froli angle is .42 rds, as 
Smee in fig.3.15 and table [. Refering to figure 3.14 there 
Meman change in the final miss distance, which is better «han 
case 0, due to a improvement in its Z componen: (see table 
i) Pres GFESul=s in @e62inai  G€G=-t5-CS miss distance of .65 
meters. 

In cas2 2, the missil2 has tha same Ac at 7=0, bué 
meeeetne cortaction for roll angie being increa 
Meee increase of target acceleration (see fig.3. 
commanded Ac reachs a larger peak value. 

Hic etc ea) value Of ene  cdumMandad roll tate 
mrgieanrs, which is larger than case 1, iue to «he increase in 
the target accelsra<ion. ii sags = 3, 13, eos e. AC 
Beeeceez ble change in the shave of che curve for 2 direc*icn, 
mm@e@mean the Y dictection the final distance is about th? same 


as in case 1. The CG-to-C3s Coweta = aspect) 25 {73 Teters. 
a 


+ 


3 
Mie larger reli tate leads to larger roll Ss seen i: 
S Ss 


fede 5.19, woere he final roll angle i 
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The effect of target accelaration can be easily seen 
in casé 3, where one can see that with the sam¢ Aco, the 
accelerations are lazrgelly increased from this point, an 
Piemeem.ssil begins its trajectory with very high roll 
Baes (see figures 3.20 and3.21). There is a change in «he Z 
component of the miss distance, thaz decreases its final 
value to .05 meters, but now the niss distance in Y direc- 
tion is made worse as shown in figur)2 3.22, which leads to 4 
final CG-to-Cs distance of 1.5 metars. In fig.3.23 one can 
se2 that the missile crosses the target with a bank angle of 
1.26 radians. 

In case 4, due to the position of target 600 meters 
under the missile, the initial commanded acceleration is 
negative and réachs the limit of -2 g's. aoe ter a) Old 
rat2 begins at a smaller value than in case 3 but increases 
@meeng the initial part 9f the flight ceaching its peak 
value at 1.75 seconds wh2n again 2S in «ae oreavious cases 
Begems «oO dscay. As the missile baaks te roil angles larger 
than 90 degrees, the accaleration goss <9 positive valués, 
Meee) in crigurss 3.24 and 3.25 Figure 3.26, shows the 


worse case among “*hese ln respect <9 tae miss distance, 


Hamby an the Z direction, and in the fiaal cgq-*o-ca 
Meseesdincs, which is equal to 4.43 nateatrs. ASO ee ana! 


m@emeangle of 3.0 radians :s the largest among all «these 


cases, aS seen in fig. 3.27. 
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Defining the projected zero 


cl 
eo 
(D 


(424) as the niss distance «he nissila 


Q 
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ta Ge wae! 


ao i ta 
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no commands. Tl can be calculatezi at et 


Miss distance 


cS 
—" 
ct 
i 
iD 


distance between target and missils pl 


ate cO The gravity or: 
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| 
fom. = z (0) + eae +2 (3234) 


which is equal in all the three f 

ieerte mn es2® > St ~tnres caSses,cne initial oo 
commanded acceleraticn is thé sSam2 aS seen in tabl2 
Considering that the control Ac nacessary “*o correc 


ewioral miss distance can be calcul2*ed as: 


BC ee ‘| rT? 2 2661 (se 35) 
Paeic 


meen 1S Cicse to the initial control Ac. 


This suggests that the initial Ac, would be chat cne 


j-- 


Meteescsalry =C Correct the initial ZEM in Z@ direction, wWaich 
agrees with reference 1. Netice aow 
M@eminitiial Ac is less than the <s 


3 A 
which agrees with the previous statemant that some 22rro0r was 
Mm 


Beemected 21n the initial part c& the computations. 
fericiie=engethe ideal iniecial “Soll aagle as dsfined 
before, ore has: 
| ZEM, 
Do See cea — == —=— (3.36) 
ideal Ze M. 


Esctecdpble ©, and Ezquzes sh 
can be verified that the nissile is b 
iE ce i Shaders: 
Beeger -han @  ., in seder to corr 
hit the target. 
Pieter G@e,in all ceases,  =2ha Missile begins its 


MeagectOry with an Aco (discussed b2ezars), and a roll rats 
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Meets Propcrtional to the targs% acceisration in they 
fees ction. As the missile rolls at decreasing roll trates, 


Pees ccemManded acceleration is changsad in ozsder for ¢t 


O 


compensate the effect of the roll angle on its 2 direction 
component.At some time when the control Pe is zéro or near 
Zero, the control Ac begins to follow a linear law, as 


suggested in eqn. 2.7: 


Ve 


Notice however that the term Aco in this equation is no more 


Ac= Aco | 1- £-| 


#he actual initial ccommanded accelsaration, but that one che 
@esstie would have if its initial roll angle was 
meses Linal . 


Such behavior defines a boundary in the cecntrol Ac 


egual to 


which is clearly seen in 3.29, where the commanded accelsra- 
t29n is bound=d by the curves of “=h2 control Ac the missils 
would have if its initial @, was 2qyual to Do ..,(OLs if no 
commanded roll was necessary to reach the target). 
Maanoudge=ne Missile commanags roll angles larger than 
=O 


omen 2d2eeal initial recll angl2s, it 1S possible 


fu 
O 
WY) 


prediction- with no no computer work- of an aproximation for 
en 
oe 


(D 


maximum acceleration the missils would experience during 


(n 


medcn, aS ctollowang (see fig. 3.32): 


ad 








ee) 
cage ip cos 
g Ordeal 
Z 6 


where Do...) cames from eqn.3.50, and Ac is a streight line 
meme Et g «3.32. 


Other interesting point is that the missile is 


Seumending tO reach roll angl2s about twice of en anes 
mie targe: is at small accelerations, and when at large 





! 
| 





acc2lerations, the missile is making a small correction on 
its roll angl2, 2S shown in table I. 

One can See from the figuztes showing miss distance, 
that the relative position of missile to target is about the 
Same in all three initial cases. Thus, the new Meal. as 
each point is <+he same. As the Coen computed at +=0 is 
smaller when the target 1S at snall acceleration 
@ymmeccion has to be larger in ordar to hit the ta 
seen in figur2 3.30. 

imeecdo we) he Che can see che LTinal miss distances, 
the final miss distance cg to cg and the final rell angle. 
Suen results show that with this iigitalized model, good 


results has been obtained. 
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Representation of the System. 


Figure 3.1 


51 





roe eT ee ee ee ee eee SEE GEE Gate’ que Gale cee eee SE ee Oc SS ee ee ee Se ee Oe ee SE ee © ee eee oe CE eee SU cee Se ge ee ee i See ee FR SS 


| 





| 


° gee SE? gg OE ggg OO ggg GE” euggy GE” guys OE eggs GE gee SEP eee GE qu ee PE que SE qqgee  qugee GSE” quay SE cee OE ge SE ge SS ge SP eee OO gee gee SS eee gee OE eee OE ee ee ee ee SP qu SS ee SE ao 


itial Roli Angle. 


3 
—=_= 


Variation of Commands with In 


Figure 3.2 


a2 





——._. ~ °° °° °° °°” 


So ee SS pce SS Go SS gy GY eee ER Ge eq ES Gee Oe ees eee ee eee eee eee eee ee” eee SE gee ee ee eee Oe ce 2 ee ee cee Ss Se em ee See SS eG eS eS 





Variation of Roll Angle - Small Angles. 


Figure 3.3 
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Variation of Roll Angle - Large Angles. 


Figure 3.4 
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Figure 3.5 Ideal Initial Roll Angles. 
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CASE 1 target acceleration -.5 g 


target inftfal pos{tfon~ CASE 1 





misstle {nitf{al position 


case 2 
target acceleration -1 g 
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Case ys dmama 
target acceleration -4 g 
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target [nftial posftfon - CASE 4 


Scenarios for Simulation. 
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Figure 3.7 


Case 0. 
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Case 0. 


Commanded Roll Rate- 
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Figure 
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Miss Distance in Z Direction- 
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IV. ANALISYS OF GAINS, SAMPLE RATE AN 


D PITCH ANGLE 


A. ANALISYS OF THE GAINS 


In chapter 3 has been developed a solution for «che 


optimal ccntrol of a system as: 


X(k*#1)= A(k) x (kK) + B(k) u (k)+ EQ 
ieewould be interesting =o check 1 the optimal gains 


reach steady state, but at the momaat char 


«t 


he sxtention for 

large roll excursions has besn introduced, and the system is 
being feed with optimal commands which are varying 2ach steno 
of time, such idea can rot be applyed. Howaver one can do 
Sem@emeensck <n th2 mod2l for small rsll excursions, which i 
acctually used to comput? the optimal commands. 

Beang this, one has the <ime history of those gains as 
in &.1 and 4.2 

Ome Can NSt1lcs in fig. 4.1 that the gain FG 
meeagee with the effect of gravity nas no ef 
Semenaed roli rats, and that FG(1,1), as shown in fig. 4.2, 
has a large Sepect. Of the tzomnmanded ac 3 
Furthermore, this gain rsachs steaiy state v 
d= Can be assumed that the gain FG wiil be equa 
steady stat2 value during all time. 


From eqn. 3.19, and assuming steady scat 


(D 


u(k)= -F(k) x(kj)-~- FG(K) g Che) 


Seas subSticuting g: 
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u(k)= -F(k) x(k)- C (4.2) 
where the second term in the right hand side is @ constant, 
and its value is exactly 2qual to *t12 value of the commended 
Meee beratlOn necessary at t=0 to correct the gravity fall of 
miemmdes2. (Or to correct the initial ZEM due t5 gravity). 

It might be supposed that one could solve the cptimal 
control preblem for the system reprasented by eqn. 4.1, Just 
considering one reduced system represented by: 


x (k+#1)= A(k) x (k)+# B(k) u (k) (4.3) 
feeereea Dias in the ccntrol as: 
u(kK)= -F(k) x(k)¢ C (4.4) 


Pumas Showed in the fdlicwing anaiysis, <tchis is not 
possibl:. 
PpemcONSceant term in <he right hand side of eqn. 4.4 is 
Sapculated as follows: 
PEeom e€qn.: 
mitcecal ZEM dus £0 gravity= ZEM = -- g <2 


C = ----- = EMg Be ee Sa ces ce ae a 
+ t? 7 
oa mg ree ccs A 


and the gains F(k) are calculated using a Ricatt+i squation 


as usual. 
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Case 5 was tested using the above consideraticns, and 
using the same scenario of case 3. 

Figures4.3 and4&.d4, show the Semis oe vO: eth > 
Some rois. On2 can see that the commanded acceleration has 
begun az same values as in case 3, but the commanded accel- 
eration reachs aop¢ak considerabl2= higher, them decreases 
does noz following a linear law, with a final at 14.7 meters 
per second syuare, being this terminal value due to the 
constan= term representing the effect of the gravity. 

Reféring to fig. 4.4, the commanded roll rate begins ac 
asame value as in case 3, but as the control Ac is too 
high, 1* reachs negative values, going *o zero almost at the 
end of the running time. This behaviour of the control leads 
to a large miss distance as seen in figure 4.5, and tabléeIt. 

Figure 4.6 shows the tim2 history of the roll angle, 


Berea ses “5 Yalues close =o 1.5 zadians. As the eccelera- 


= 


memeat ~hiS point is larger than the correct value, the «the 
corrections are excessive and the roll angledecreases at «he 
emaeeont she running time to the valuz of .12 radians 
mods, One® can see that the gain due «+ 2 ew gaat! S 


O° 
S@eeesestation can net be replacsi by its st 


(D 


ady stats¢ 
a 


@Ome 12 


” 
a 
O 
+ 
Ce: 


meme. rhs kind of simplification san thu 
S 


Be EFFECT OF THE SAMPLE RATE 


Mimeughcut all the simulati:os 2 sampl= osricd cf .05 
S=conds has been used.In this section 4&4 brie=i study of che 
pmeeone Of the change cf this samples rate is given. 

Two best cases have been ssléecced ¢ illustrate whe 
effect of the sample period. 

The first case, case 6, has been run with a sample 


meerod Of 41 SécondsS and corsits > the same scenario as 
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As one can see in figure4.7 andj taolsIIii, theme 1S) 10 
noticable change in the commanded acceleration, but the 
commanded rcll rate begins at small2r value than in case 3, 
as is seen in fig.4.8 This initial decreasé in roll rate, 
leads to a large miss distanca ina ¥ direction as shown in 
fig.4.9, and to a small v2alue of roll angle (see fig.4.10). 

The second cas2, cas2 7,has 32 smple period cf 0.025 
seconds. There is no noticeable change in the tin 
@weene control Ac as shown in fig¢g.4.11. The comman 
rat? begins at a higher value than in case 3 as Sh 
Mmeerwntch Leads to a final miss distance in Y direction of 
2.22 metsrs andin Z direction smaller than case 6 (see 
fig4.13). Figure4.14 shows that the final roll angle is 
increased and the missile cross the target with 1.28 radians 
and with a CG-to-CG distance of 2.5 meters. 

In both cases, the miss distaat2 was increased over «he 
nominal value obtained with a sampling rate of 0.05 seconds. 
Thus, i= would appear that there is an optimal value for the 
Sampling ra+t?, which may be conected wit the geomezry of 


on 
the scenario and with time *0 go. 


C. EFFECT OF THE INITIAL PITCH ANGLE 


[mens 2mMmportent at this step > remenber that througout 
this work has we have peen discussia: a dimensional model, 
Mmemescthere is no infcrmetion on th= X coordinate, so it is 


impossible <=c verify the behavior o 


th 
ct 
= ie 
(Dp 

re) 
4 
ct 
Q 
etn 
by 
ig 

Q 
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iD 
a 


iments work, Sincerity alleen 


(v 
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Db 
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4 
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fs 
Ut 
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Q 
WwW 
“J 
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O 
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ct 
oe 
(D 


initial angle 9 was equal to zero, this valus has been kept 
MemecOnStant dtring all time, and considering that whithout 
@avereniLOrmation of a third dimension 1t was not possible to 
Senrrece the tine to intercept, this time was also kept 


h 
constanz and equal tc the nominal value of 5 seconds. 


Se 





ioeace chat <chis assumptzon is likely ¢ 


O 
Smee nas the horizonzal initial distance fron targs <> 
missile compared with the initial vertical distance between 
target and missile large enough in Ogicu, <O have 
small angles. 
The question that rises is, howecoOGrdecnss Paech engle 
affect the system if it was net smali? 


As seen in fig.4.15, the missile velocity in the xX 
direction would be: 


te ¥,, cos8 + g cos sind t (4.5) 


Memeerr nes an effect not only from the pitch angle, repze- 
Samered by the cos , but also from the gravity's accelera- 
meee wniich will ie@ads to a different tine to intercept. 
Considering «he same physical sctinaris as in case 4, bu 
Seeeaing the initial pitch angle, iewoOmes= “tO have the 
Meese pointing t> the target (s22 f2g.4.16), and keeping 
eeeommassiles velocity cf 1000 m/sec in the X direction, ons 
has a ccempletly different geometry of the problem as seen 


from the flight path reference frame. 


With this new situation (see fiz. 4.16), case 8 has been 
iain < Figures 4.17 and 4.18 show ‘that now the missile is 
te jee ee as 


Setanding lazrj= positive accelerations, and the cr 
meempeg. ning of the flight is too high, g 

very small pariod of time. The 
#3g.4.19 are increased in the initial p 
@eectne missile correc 

Sats te ee gaol 


7 
reach a final cCG-to-CG miss dista 
ValcmmesS OSCCelsc Oty in the 


aegeee, due £5 he large cont= 
begining and becames constant with avaiue of .57 radians 
(see f2g.4.20 and tabl2 IV). 
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Notice that the high values of acceleration needed are 


in som? part due t0 a vertical component of targexts 
VemgGeniey, which is seen From the flight path frame as the 
target was manouvering in the 2 PESece2On Ween Cons tan~ 


Welbocity. These large accalerations leads to roll zréetes too 
Baege for che physical integrity cf the missile. Thi 
em@a@enat hough althoug the gcod results obtained, if 


with case 4, they are not practical. 


Mm 


Mpeorde=> te get rid sf the vertical manousver of th 


target, casé 9 has been run. In this case, zhe Scenario is 


ia 


jmemoame aS bafors with the <argez also pointing down, wit 
the same pitch angle as the missil=, and has a X velccity 
equel *o the previous cass (s2#2 fig.e4.16). 


Peom fig.4.21, one can see that tha decrttzase in the 
control Ac is substantial if compar: i-escase. CG. hus she 


u 
Semmended roll rate is still toc larg? as shows fig.4.22 The 
feewema.StOry Of the miss jiistance is batter, resulting 
final cCG-to-CG distance of abcur 50% £ c2s¢e 8 =(se23 
foment) 23). Winon mores ss )g ier = Ss aos Oscttatory .asS Sse 
meme. 24 and the missiles cross the targez with 1.3 radians 
in roll. (see tablerv). 


Mae results cbtained in the two previous ca 


(i 
iD 
"A 
= 
{n 
(> 
Q 
OQ 
Mm 
72) 
cf 
W 


Smeethe algorithm developed in this work cculd be readly 
eeeested tO alr-to-surface missiles. fe whew Laicct, . the 
Se=nNatic would be favorabl2= to ths missile than in 2ither 
Peemeous cases, Since on2@ can csorsides thaz the *arge=z cou 
be essencially stationary in comparison with the miss 
spe2d. 

Case 10 has been run with this en eae snd the scen- 

sf 


h 
Seio as in fig. 4.16. Ths target i 
b 


acceleration, and the missile be Pleghe 600 met er 

meeve i+, with the same initial pitch angle as before. 
Pgete25 SNOVvS =he tins History or he commanded accel- 

eration, where one can sé2 that as there is no roll rate to 
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be commanded, the acceleration is following a straigh= lins 
with very resoanable values. The miss distance is show in 


figure4.27, which shows the final CG-to-CG distance of .31 


meters. 
Based in the results of these tests, one can ses that 
there will be some effect of the pitch angie on the miss 


Meerenc>, not only due to its effect on the time to inter- 
cept, but also because at the moment that thera is a pitch 
angle differeat from zero, even if the targerz is 
meeame Level, in the tlight path frame &@ component of the 
target's velocity will show up leaiing the miss 
command large accelerations and roll rates. AiehOy Gg it has 


results harm the performance in an air-to-air missiles, in 


tha case of air-to-surface missilas, when the target has 
been considered with no motion, good results has been 
obtained. 


feeorr ECT OF TIME TO INTERCEPT 


Mmemeene Si .mulation of case 1 and 2 in chaptar 3,1* has 
been observed that when the target was az small accelera- 
mons, <=he missile did larger corrections on its roll angle, 
meee respect «=o the ideal initial roll angle, than wher the 
Beeget Was With large accelerations. One can think that aust 
be scme kin ef compromise batween Vv 
targe+ and missile (WHICH weet Gor 
ieeercep*), and the relative position be 
will affect the miss distance. 

mimorde= £0 dO a brief analzsyS on this, case 11 and 12 
has béeén run. 

Pyecase 11, the scenario of case 2 has been kept, with 
Biee excep ~ion that the missile's velocity was change +o 2000 


Mesec., Which means that Ti was changed to 2.5 seconds. 
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Figure 4.29 shows that the acceleration is largely 
iicceased due to the small time required to correct *h2 ZEM 
and the commanded rcll rate is almosst twice of case 2 (s 
fgeeed.30). Fhe final mists distance is more th2n four tin 
the value obtained in case 2, as s22n in fig. 4.31 and tab 

n 


V. The final roll angle is about half as in cés@2@ 2, si 


Q aA F- 
iD 


the projected final position of tn= target in the Y direc- 
tion is less than in case 2 (see fiz. 4.32). 

Case 12 was cun with a scenaris less favorable to the 
missile, where all the conditions of the previous cas¢ was 


kept, except the target position that has been increased ¢ 


O 


200 meters above the missile. 

Now, the sommanded accaleration ars much larger, with a 
initial Aco of 103 m/séec2, being almost impossible to see 
the difference of the *ime histor Eeemereceslsratzon from 
@mematcraigh= line, as shown in fig. 4.33 The commanded roll 
Tate is small, about che same aS in case 2 (see Fig. 4.34). 
Piemenhange in the miss distance is noticeable, with a final 
Wemeg~cG distance of 4.8 meters as in table V, and figure 
4.35. The final roll angl2 explain tha shape of the acceler- 
eel Curve, Saimeecewlen the small @rotl anele as Shown in 
mage +. 30, the system is behaving 2s for small roll éexcur- 
2aonS. 

Meeece that from this analyzis, one has t> realizes «hast 
there is seme kind of envelove a she 2-D syste is 
valid. And in ordér te datermine chis envelope, one has to 
take in acccunt not ecnly the geomezrry of «he 
Beso tn] time to intercept, which is determin 
mien celation oF velccities of missile and <ar 


Beech aigle *50. 
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Figure 4.1 
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Figure 4.4 
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TABLE II 
Results from Biased Control 


miss 

miss 

distance distance 

Y direction ¢ Z direction 
(m) {m) 
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Figure 4.21 
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The scope of the present work was the development of an 
+ 


W 
Seeemal digital control Ce iewearoliogsoume 4 Dank=<O turn 
missils. 

A two dimensional modal, as sug2st2d in reference 1, was 
meepeed. Atcver the digitalization of the comtinuous model it 


-_~ 


4 


WwaS necessary *0 solve a modified Ricatti equation since 


~ 


zt 

the state equation there was a third term representing «hs 
gravity's effect. The approach that has been adopted is new, 
and although good results were obtained for *h2 scenarios 
considered in this work, is necessary that she algorithm bs 
further tested and evaluated in Similar problems du 
Mee 1 cy. 

The optimal was solvei with aa a 
small angles. This condition was later rela 
Bemel angles could be analyzed. 


Mme iS dirticult to compare the praesent work wit révious 


v 
eB 
te 
{nh 
ct 
ay 
(Pp 
1 
re 


results since Stallard has indicated his orig- 
aial paper, and further works in this area was net found. 
However some comparison with Stallard work is possible 
The commanded acceleration cf th2 nissils ar2 such as to 
eorrect the ZEM at each point, this agrees with tha« refer- 
ence.There is a proportional relationship between the 


am 

commanded recll rate and *he commandsd acceleration , and the 
Bemmandsd rcll rate is proportional to the defined Cee at 
each point, which again agrees with 

The algorithm developei in this work requirss extensive 
computaczion at each step, and is slear tha«+ some softwares 
O9t2Miza ion will be needed.The motivation for considering 
the constant pesos e> Ss oa tee ailemm co Gravity Was <0 
decrease this ccmputational burden. This apprsach however 


resulted in unacceptable miss distanc:3. 
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Mio Me ameocm or Of Invest ~¢dat1On that conld reflect on th: 
period available to the computer to perform its calculations 
was a change in +he sampl2= rate. Iwo different sample rates 
were investigated, both lead to larger errors than the 
nominal period of .05 seconds. A detailed study on this 
issue is ‘left as sugestion for future works, Since some 
optimal value of the sampling rate is clearly indi 

It is important to kasp in mind that the mo 


is two dimentinal, while the actual problem is ¢ 
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te! 
Do 
O 
4 
fu 
Mm 
ry 
ct 
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Slonal, thus some brief studies we2rs conduc 
check the region of validity of the 2-D. 

In the analysis of the pitch angle, one can see *hat i 
necessary +o have small variations in pitch in corde 
approximate it as a const2nt. However, at the mon 
Somes angle is different from zero, as explained in 
4, it 1s possible to have in the flight oath rs 


c 

a target manouver in the Z direction that will lead «o lards 
S 
3 


acceleration tommands, leading tn2 missile ¢ large miss 
distances, when considéring a movable target. When the 


resent syStem was tested tested against fixed targets, the 


b 


mol tS were guite gcod, this sugjy2sts the applicaticn of 
*ha2 model in air-to-surface missiles. 
Further investigation were mad> on th 
ao= GO. AS expected, decreased time? to go, results ia 
increased miss distance. A detailsd a 
Somplex Scenarios is needed in order to 
memect Of tims to go. 
MEsGummc ss WOULG be intereszcing <5 esxcend the model to 
Peies]e¢ dimensions and include <ne effects of lLags on the 
system in future works. Beta, ih appendix A, the 
computer model used in this work is :nclosed. Some improve- 
méneS in this program can be done, mainiy in the data intro- 


mies. 0ON, and 12 somes Optimization of the running time. 
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APPENDIX A 
FORTRAN PROGRAM 


These appendix provides a listing of the computer 
program used in the present study. 

Since the routines used ar2 non-IMSL, and a smail change 
+o double-precision was necessary, they are also being 


provided. 
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